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ABSTRACT 


An  adaptive  focusing  microwave  planar  phased-array  hyperthermia  system 
for  improved  heating  of  superficial  tumors  is  experimentally  investigated.  A  com¬ 
mercial  microwave  hyperthermia  planar  phased-array  antenna  system  at  the  Uni¬ 
versity  of  California  at  San  Francisco  (UCSF)  Radiation  Oncology  Department 
has  been  modified  to  implement  a  gradient-search  adaptive  focusing  algorithm. 
Adaptive  focusing  measurements  in  a  muscle-equivalent  liquid  phantom  with  a 
915-MHz  microwave  planar  phased-array  hyperthermia  system  comprising  16  in¬ 
dependent  amplitude/phase-controlled  waveguide  antenna  elements  (Microtherm- 
1000,  Labthermics  Technologies,  Inc.,  Champaign,  IL)  are  reported.  An  electric- 
(E-)field  feedback  probe  detector,  fabricated  at  UCSF,  measured  the  E-field  gener¬ 
ated  by  the  hyperthermia  phased  array.  A  method  of  steepest-ascent  gradient-search 
feedback  algorithm,  implemented  in  software,  controls  the  hyperthermia  array  phase 
shifters  and  focuses  the  transmitted  radiation  beam.  In  10  to  15  iterations,  the  mea¬ 
sured  phantom  data  indicate  rapid  convergence  of  the  adaptive  focusing  algorithm 
and  significant  increase  of  the  focal  region  field  strength  due  to  the  adaptive  fo¬ 
cusing.  Two-dimensional  El-field  radiation  pattern  measurements  were  collected  by 
scanning  the  El-field  dipole  probe  antenna  inside  the  muscle- equivalent  liquid  phan¬ 
tom.  The  measured,  adaptively  phase-focused  radiation  pattern  data  indicate  a 
maximum  useful  heating  depth  of  about  3  cm  in  a  muscle-equivalent  sugar/saline 
phantom  having  dielectric  losses  of  3  dB/cm. 
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1.  INTRODUCTION 


Adaptive  array  antennas  are  well  known  for  their  ability  to  improve,  in  real  time,  the  perfor¬ 
mance  of  communications  and  radar  systems  [1-4].  Recently,  adaptive  array  techniques  have  been 
applied  to  the  medical  use  of  radio  frequency  (RF)  hyperthermia  for  deep-seated  tumor  therapy  [5- 
12].  With  an  adaptive  RF  or  microwave  hyperthermia  array,  it  is  possible  to  automatically  control 
the  electric-  (E-)field  at  multiple  positions  within  the  target  body  [5-12],  The  E- field  radiated  by 
a  hyperthermia  phased  array  can  be  minimized  (nulled)  and  maximized  (focused)  at  desired  target 
positions  by  adaptively  adjusting  the  transmit  amplifiers  and  phase  shifters  of  the  hyperthermia  ap¬ 
paratus.  Multiple  adaptive  E-field  nulls  and  adaptive  focusing  were  experimentally  demonstrated 
at  the  State  University  of  New  York  (SUNY)  Health  Science  Center  with  a  modified  commer¬ 
cial  RF  hyperthermia  phased-array  system  and  phantom  target  measurements  [9-12].  An  annular 
phased-array  antenna  (BSD-Medical  Corp.,  Sigma-60  [13-14]),  operating  at  «  100  MHz,  was  used 
[9-12].  This  report  investigates  a  915-MHz,  planar  phased-array  antenna  hyperthermia  applicator 
via  adaptive  focusing  experiments  at  the  University  of  California  at  San  Francisco  (UCSF). 

Treating  a  patient’s  malignant  tumors  is  often  a  difficult  task.  The  objective  of  the  treatment 
is  to  reduce  or  completely  remove  the  tumor  mass  by  one  or  more  modalities  available  at  the 
treatment  facility,  commonly  surgery,  chemo-,  and  x-ray  therapy  [15].  One  particular  method  used 
by  itself  or  in  conjunction  with  another  is  “tissue  heating”  or  hyperthermia  [15-19],  a  form  of  high 
fever.  A  controlled  thermal  dose  distribution  is  required  for  hyperthermia  to  be  a  valuable  therapy. 
Typical  localized  hyperthermia  temperatures  required  for  therapeutic  treatment  of  cancer  are  in  the 
42.5°  to  45°C  range.  Normal  tissue  should  be  kept  at  temperatures  below  42.5°C  during  treatment. 
The  most  difficult  aspect  of  implementing  hyperthermia  with  either  RF  or  acoustic  (ultrasound) 
waves  is  producing  sufficient  heating  at  depth.  Multiple-applicator  RF  hyperthermia  arrays  are 
commonly  used  to  provide  a  focused  main  beam  at  the  tumor  position.  A  focal  region  should  bp 
concentrated  at  the  tumor  with  minimal  energy  delivered  to  surrounding  normal  tissue.  As  the 
hyperthermia  antenna  beamwidth  is  proportional  to  the  wavelength,  a  small  focal  region  suggests 
that  the  RF  wavelength  be  as  small  as  possible.  Due  to  propagation  losses  in  tissue,  however, 
RF  penetration  depth  decreases  with  increasing  transmit  frequency.  Typically,  noninvasive  phased 
arrays  operating  close  to  ICO  MHz  are  suggested  for  deep  seated  and  915-MHz  arrays  for  shallow 
tumor  heating.  One  of  the  major  problems  in  heating  a  tumor  with  a  hyperthermia  antenna  is 
the  formation  of  undesired  hot  spots  in  surrounding  tissue,  which  often  produces  pain,  burns,  and 
blistering  in  the  patient,  requiring  terminating  the  treatment.  Thus  techniques  for  reducing  hot 
spots  are  necessary. 

A  typical  clinical  RF  or  microwave  hyperthermia  treatment  consists  of  several  two-hour  ses¬ 
sions  spread  over  a  period  of  weeks.  The  patient  is  prepared  during  the  first  hour  of  the  session — 
temperature  and  vital  signs  monitoring  instrumentation  axe  interfaced  with  the  patient,  the  hy¬ 
perthermia  antenna  microwave  signal  generator  is  turned  on,  and  the  tumor  is  heated  by  radiated 
electromagnetic  energy.  The  tumor  should  be  heated  for  «  45  min  at  a  temperature  of  42.5°  to 
45°C.  Prior  to  the  actual  delivery  of  microwave  energy  to  the  tumor,  «  15  min  (maximum)  is 
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available  to  manually  or  adaptively  shape  the  microwave  E-field  distribution  within  patient  target 
tissue.  Current  clinical  operation  of  commercial  hyperthermia  phased  arrays  allows  limited  man¬ 
ual  control  of  the  array  transmit  element  amplitude  and  phase.  Some  improvement  in  the  E-field 
distribution  can  be  achieved  by  this  manual  trial-and-error  method,  but  automatic  adjustment 
techniques,  offered  by  adaptive  arrays,  are  desirable  and  possibly  offer  better  E-field  distributions. 

Several  journals  have  published  special  issues  on  acoustic  and  electromagnetic  hyperthermia 
treatment  of  cancer  [20-22].  Studies  have  been  conducted  to  produce  improved  therapeutic  field 
distributions  with  hyperthermia  phased  arrays  [23-42].  Phase  control  can  synthesize  improved  RF 
radiation  patterns  without  adaptive  control  of  the  transmit-array  weights  [23-28].  Array  transmit 
weights  (defined  here  as  amplitude  and  phase  states  associated  with  a  transmit  channel)  can  be 
adaptively  controlled  to  maximize  tumor  temperature  (or  microwave  power  delivered  to  the  tumor), 
while  minimizing  the  surrounding  tissue  temperature  (or  RF  power  delivered  to  the  surrounding 
tissue)  [29-42],  Ultrasound  receive  feedback  sensors  [41,42]  and  a  pseudoinverse  pattern  synthesis 
method  [43]  generate  multiple  focal  points  for  a  phased-array  hyperthermia  applicator.  Previous 
reports  [5,9]  and  several  conference  articles  [6-8]  investigate  the  theoretical  benefit  of  using  adaptive 
nulling  with  noninvasive  auxiliary  dipole  sensors  to  reduce  the  E-field  intensity  at  selected  positions 
in  the  target  body  while  maintaining  a  desired  focus  at  a  tumor.  Multiple  adaptive  E-field  nulls 
were  used  to  show  a  theoretical  reduction  in  hot  spots  for  a  homogeneous  elliptical  phantom  target 
surrounded  by  a  water  bolus  and  hyperthermia  ring  array  [6-8]. 

Microwave  antenna  applicators  are  often  considered  for  inducing  localized  hyperthermia  to 
superficial  tumors  located  on  the  chest  wall  and  head  and  neck  regions.  The  power  deposition 
patterns  (also  referred  to  as  the  specific  absorption  rate  or  SAR)1  of  single  aperture  applicators 
are  often  limited  to  effective  heating  well  within  the  aperture  boundaries  [44-47].  With  915-MHz 
microwave  radiation,  waveguide  applicators  have  a  limited  heating  depth  to  less  than  3  cm  and 
lateral  heating  dimensions  of  3  to  5  cm,  which  is  smaller  than  needed  to  cover  many  tumors 
[48].  Heating  depth  is  limited  to  only  about  3  cm  because  at  915  MHz  the  dielectric  loss  in 
tissue  is  typically  3  dB/cm.  A  predictor  of  good  local  heating  control  is  to  use  the  50%  iso-SAR 
coverage  (laterally  and  with  depth)  of  the  tumor  [49].  Several  methods  for  improving  the  SAR 
heating  patterns  of  single-aperture  hyperthermia  devices  are  under  investigation.  For  example, 
microwave  absorbing  saline-filled  boluses  can  shape  the  heating  pattern  of  a  waveguide  applicator 
as  demonstrated  in  Sherar  et  al.  [50].  Mechanically  scanned  microwave  applicators  have  been 
used  to  shape  the  power  deposition  pattern  [48,51,52].  Phased-array  applicators  can  shape  the 
power  deposition  pattern  produced  by  planar  [48,51,53,54]  and  conformal  arrays  [55,56].  Studies 
to  improve  penetration  depth  with  phased  arrays  by  controlling  the  phase  and  amplitude  of  each 
array  element  have  been  conducted  [53,54,57]. 


!SAR  =  \o\E\2 / p,  where  E  is  the  E-field,  a  is  the  electric  conductivity  of  the  tissue,  and  p  is  its 
density. 
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The  Microtherm- 1000  (Labthermics  Technologies,  Inc.,  Champaign,  IL)  hyperthermia  mi¬ 
crowave  planar  phased-array  antenna  system  used  in  these  experiments  is  shown  in  Figure  1.  By 
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Figure  1.  The  915-MHz  microwave  planar  phased-array  hyperthermia  system  used  m 
experiments  i‘.  the  UCSF  Hyperthermia  Laboratory. 


transmitting  close  to  the  patient  with  a  planar  phased  array,  it  is  possible  to  obtain  superficial 
heating  over  large  («  15  x  15  cm)  as  well  as  small  («  3.5  x  3.5  cm)  areas,  depending  on  the  array 
amplitude  illumination  function.  In  Figure  2,  the  Microtherm- 1000  uses  a  15-  x  15-cm  planar  ar¬ 
ray  with  16  square  waveguide  elements  operating  at  915  MHz.  The  16  independent  variable  power 
amplifiers  drive  the  waveguides,  and  each  of  the  16  active  channels  has  an  electronically  controlled 
variable- phase  shifter  to  focus  the  array.  A  cool- water  bolus  between  the  patient  and  the  phased 
array  prevents  excess  heating  of  the  skin  surface.  The  bolus  is  filled  with  circulating  deionized 
water,  which  has  a  very  low  microwave  propagation  loss  of  about  0.3  dB/cm. 

A  candidate  adaptive  focusing  algorithm  proposed  by  MIT  Lincoln  Laboratory  for  use  in 
hyperthermia  is  a  gradient  search  based  on  maximizing  the  signal  received  by  an  E-field  sensor 
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Figure  2.  The  915-MHz  microwave  planar  phased-array  applicator  used  m  the  Labiher- 
mics  Microtherm- 1000  system. 


positioned  within  a  tumor.  UCSF  has  recently  evaluated  (with  experiments  in  phantoms  and 
animals)  a  systematic  search,  iterative  focusing  algorithm  for  phased-array  hyperthermia  [58].  Note 
that  the  E- field  probe  used  in  the  UCSF  measurements  is  a  short  dipole  with  a  semiconductor  diode 
detector  [59-62].  A  particular  concept  of  adaptive  focusing  as  it  applies  to  hyperthermia  follows. 

Recently,  a  technique  called  focused  near-field  adaptive  nulling  [63-70]  was  investigated  to 
test  the  performance  of  radar  and  communications  systems.  A  calibration  probe  antenna  focuses 
the  main  beam  at  approximately  one  diameter  D  of  the  phased  array.  Additional  (auxiliary)  probe 
antennas  are  located  in  the  sidelobe  region  of  the  phased-array  quiescent  radiation  pattern.  Adap¬ 
tive  nulls  are  formed  in  the  direction  of  these  auxiliary  probes.  Computer  simulations  show  that 
a  radar  system  can  be  tested  by  using  the  focused  near-field  nulling  technique  [68,70];  experimen¬ 
tal  measurements  of  focused  near-field  adaptive  nulling  have  been  performed  [69].  For  the  radar 
system,  interference  signals  entering  through  uncontrolled  sidelobes  can  reduce  the  signal-to-noise 
ratio  (SNR);  adaptive  nulling  [1-4]  counteracts  this  possible  degradation.  In  the  hyperthermia 
system,  high-transmit  antenna  E-field  strengths  can  create  hot  spots  in  the  target  tissue.  These 
high-temperature  regions  can  possibly  be  alleviated  with  adaptive  nulling  of  the  E-field  at  specified 
positions  on  the  surface  of  or  inside  the  target.. 
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Consider  a  potential  clinical  application  of  the  proposed  adaptive  hypertheimia  system.  The 
Microtherm- 1000  currently  does  not  supply  El- field  probes  to  monitor  clinical  hyperthermia  treat¬ 
ments;  however,  as  current  measurements  indicate,  an  E-field  probe  could  be  added  to  this  system 
to  provide  feedback  signals  to  the  adaptive  algorithm.  In  theory,  with  16  independent  transmit 
channels,  the  El- field  radiation  pattern  can  be  controlled  at  up  to  16  points  (in  the  radiation  field), 
using  a  feedback  signal  measured  at  each  point. 

The  remainder  of  this  report  is  organized  as  follows.  Section  2  discusses  the  concept  for 
an  adaptive  focusing  hyperthermia  system  and  describes  in  detail  the  gradient-search  feedback 
algorithm  implementing  the  adaptive  array.  Section  3  describes  the  materials  and  methods  used 
to  obtain  the  measured  adaptive  focusing  results  in  a  phantom  target.  Section  4  gives  measured 
results  for  an  adaptive  focusing,  hyperthermia  planar  phased  array  operating  at  a  CW  frequency 
of  915  MHz.  The  measured  received  microwave  power  distributions  at  the  intended  focus,  before 
and  after  adaptive  focusing  with  a  short-dipole  field  probe,  are  presented.  Limited  two-dimensional 
El- field  radiation  patterns  are  also  shown.  Section  5  presents  conclusions. 
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2.  THEORY 


2.1  Adaptive  Focusing/Nulling  Hyperthermia  System  Concept 

The  concept  of  an  adaptive  focusing/nulling  hyperthermia  system  is  shown  in  Figure  3.  The¬ 
oretically,  to  generate  the  desired  field  distribution  in  a  clinical  adaptive  hyperthermia  system, 
receiving  sensors  are  positioned  as  closely  as  possible  to  the  focus  (tumor  site)  and  where  high 
temperatures  are  to  be  avoided  (such  as  near  vital  organs  and  scar  tissue).  A  noninvasive  adaptive 
nulling  system  is  achieved  by  placing  auxiliary  sensors  1,2,”  •  ,Naux  on  the  target  skin  as  shown. 
(With  915-MHz  microwave  radiation  as  considered  in  this  report,  hot  spots,  if  present,  are  expected 
to  occur  near  the  target  surface.)  The  null  zones  centered  at  each  auxiliary  probe  naturally  extend 
into  the  target  region  to  eliminate  undesired  hot  spots.  The  width  of  each  null  zone  is  directly 
related  to  the  strength  of  each  null  (sometimes  referred  to  as  the  “amount  of  cancellation”),  which 
is  directly  related  to  the  SNR  at  the  sensor  position.  A  low  SNR  produces  a  small  amount  of 
nulling,  a  high  SNR  produces  a  large  amount.  Resolution  (minimum  spacing  between  the  focus 
and  null  positions)  is  normally  equal  to  the  half-pcwer  beamwidth  of  the  antenna  and  is  enhanced 
somewhat  by  using  weak  nulls  whenever  the  separation  between  the  null  and  focus  is  closer  than 
the  half-power  beamwidth. 

Initially,  the  hyperthermia  array  is  focused  to  produce  the  required  field  intensity  at  the 
tumor.  An  invasive  probe  is  required  to  achieve  the  optimum  focus  at  depth.  To  avoid  hot  spots,  it 
is  necessary  to  minimize  the  power  received  at  the  desired  null  positions  and  to  constrain  the  array 
weights  to  deliver  a  required  amount  of  transmitted  or  focal  region  power.  The  adaptive-array 
weights  (with  amplitude  A  and  phase  $)  are  controlled  by  the  gradient-search  algorithm  to  rapidly 
form  the  focused  beam  before  a  significant  amount  of  target  heating  takes  place.  With  this  adaptive 
technique,  it  should  be  possible  to  avoid  hot  spots  and  maintain  a  therapeutic  thermal  dose  at  the 
tumor. 

2.2  Adaptive  Transmit- Array  Formulation 

Consider  a  hyperthermia  array  with  N  antenna  elements,  where  N  =  16  for  the  Microtherm- 
1000  system.  The  input  signal  to  each  of  the  N  array  elements  is  obtained  from  the  amplitude  and 
phase-adjusted  signal  distributed  by  a  power  divider  network.  The  number  of  adaptive  channels  is 
assumed  to  be  equal  to  the  number  of  transmit  antennas  A.  Lettu  =  (tui,u;2,  ••  ',w  yv  )T  (superscript 
T  =  transpose)  denote  the  adaptive  channel  weight  vector  as  shown  in  Figure  3. 

Commonly,  the  weight  vector  is  constrained  to  deliver  a  required  amount  of  power  to  the 
hyperthermia  array  or  the  tumor.  For  simplicity  in  the  experimental  adaptive  hyperthermia  array 
control  software,  the  weights  are  constrained  such  that 

/v 

X>n|2  =  A-  ,  (1) 

n=l 
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Figure  3.  Adaptive  focusing/nulling  hyperthermia  system  concept. 
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where  jwn|  is  the  transmit- weight  amplitude  for  the  nth  adaptive  channel  and  K  is  a  constant.  To 
generate  an  adaptive  phase  focus,  a  gradient-search  algorithm  controls  the  transmit  weights  (phase 
shifters). 

2.2.1  Gradient-Search  Algorithm 

Gradient-search  algorithms  are  commonly  used  in  adaptive- array  applications,  particularly 
where  the  channel-to-channel  correlation  cannot  be  calculated  or  measured.  With  a  gradient  search, 
only  the  output  power  of  the  receiver  channels  is  measured  and  used  as  a  feedback  signal  to  the 
algorithm.  A  wide  variety  of  gradient  searches  exists  [71-76];  this  report  uses  a  standard  method 
of  steepest  descent  [71,77]. 

In  the  Microtherm- 1000  system  at  UCSF,  only  the  received  E-field  power  is  measured,  hence 
it  is  appropriate  to  consider  a  gradient-search  algorithm  to  control  the  E-field  power  at  selected 
positions.  Gradient  search  controls  the  transmit  weights  iteratively  so  that  the  microwave  signal 
received  by  the  probe-array  is  maximized  (focused).  Transmit-array  phase  shifters  are  adaptively 
changed  in  small  increments  and  probe-array  output  power  is  monitored  to  determine  phase  settings 
that  increase  the  output  power  most  rapidly  to  a  maximum.  The  mathematical  formulation  for  the 
gradient  search  is  straightforward  [71,77]  and  is  next  described  in  the  context  of  hyperthermia  [9]. 
Although  the  formulation  is  given  as  a  maximization  problem,  the  equations  are  readily  converted  to 
the  minimization  problem.  Considered  here  is  the  general  mathematical  formulation  of  amplitude 
and  phase  dithering;  however,  only  the  phase-dithering  approach  is  applied  in  the  adaptive  focusing 
experiments. 

The  summation  of  power  received  at  the  E- field  probe(s)  is  denoted  prec.  The  adaptive- array 
signal  enhancement  ratio  C  is  defined  here  as  the  ratio  of  the  summation  of  probe-received  power 
after  adaption  p„  to  the  summation  of  probe-received  power  before  adaption  p^;  that  is, 

C=^  .  (2) 

Pb 

Consider  now  J  sets  (or  iterations)  of  JV  transmit  weights  that  are  applied  to  an  adaptive 
hyperthermia. phased-array  antenna.  In  terms  of  adaptive  focusing,  the  optimum  transmit- weight 
settings  (from  the  collection  of  J  sets  of  N  transmit  weights)  occur  when  the  total  power  received 
by  the  auxiliary  probe  array,  denoted  prec,  is  maximized.  For  notational  convenience  let  a  figure  of 
merit  F  denote  prec  and  employ  a  method  of  steepest-ascent  gradient  search  to  find  the  optimum 
transmit  weights  to  maximize  F;  that  is, 

Fopt  =  max(Fj)  j  =  1,2,* ••, J  .  (3) 

Ideally,  in  the  adaptive  focusing  gradient  search,  power  should  increase  after  each  iteration, 
that  is,  F2  >  F\,  F3  >  F2,  F4  >  F3,  •••,  Fj  >  Fj-\.  Due  to  system  noise,  however,  this  ideal 
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behavior  may  not  always  be  achieved.  Assume  that  there  are  N  complex  transmit  weights  in  the 
hyperthermia  phased  array  as  suggested  by  the  amplitude  and  phase  scatter  diagram  depicted  in 
Figure  4.  The  nth  transmit  weight  in  the  jth  configuration  (or  iteration)  of  transmit  weights  is 
denoted 

t vnj  =  Anje3*n>  ,  (4) 

where  Anj  is  the  transmit- weight  amplitude  distributed  over  the  range  Am,„  to  Amax  and  is  the 
transmit-weight  phase  distributed  over  the  range  $min  to  $mai.  The  goal  is  to  find  the  values  of 
amplitude  and  phase  for  each  of  the  N  transmit  weights  such  that  F  is  maximized,  and  an  adaptive 
radiation  pattern  peak(s)  is  formed  at  the  auxiliary  sensor  position(s). 


1(7121-2 

<J> 


Figure  4-  Scatter  diagram  in  rectangular  coordinates  for  the  amplitude  and  phase  of  the 
transmit  weights  in  a  hyperthermia  phascd-amy  antenna.  The  nth  transmit  weight  in  the 
jth  configuration  of  transmit  weights  is  denoted  w„j. 
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Assuming  an  initial  setting  of  the  N  transmit  weights,  the  weights  are  adjusted  by  dithering 
them  until  the  optimum  F  is  achieved.  The  goal  is  to  find  the  collective  search  directions  for  the 
N  transmit  weights  such  that  F  increases  most  rapidly;  that  is,  transmit  weights  are  selected  so 
that  the  directional  derivative  is  maximized  at  (A j,$j),  where  A:  and  are  the  amplitude  and 
phase  column  vectors,  respectively. 

The  directional  derivative  of  Fj  is  expressed  in  terms  of  the  amplitude  and  phase  changes  of 
the  transmit  weights  as 


(5) 


where  d  means  partial  derivative  and  r^n;,  r$nj  are  the  (A,$)  directions  for  which  Fj  is  increasing 
most  rapidly.  The  directions  rAnj,  r$nj  are  constrained  by 


N 


=  1 

n=l 


(6) 


The  goal  is  to  maximize  D(Fj)  subject  to  Equation  (6).  Using  Lagrange  multipliers,  construct  the 
Lagrangian  function 

N  ac\  a  v.  N 

Lj  =  5Z( pT^nj  +  a^^nj)  +  G[1  -  +  rjn>)]  ,  (7) 

n=l  OAn>  09 n=l 

where  G  is  a  constant  to  be  determined.  The  requirement  that  Lj  be  an  extremum  implies 


or 


dLj 

drAnj 


dLj 

dr$nj 


dAn} 


d*nj 


2GrAnj  =  0  ,n  =  1,2,  *  -  - ,  TV 

(8) 

2Gr$„j  =  0  ,n=l,2 , 

(9) 

1  dfj 
TAnj  2 G  dAnj 

i  dFj 

r*nj  "  2 G 


(10) 

(11) 
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Squaring  Equations  (10)  and  (11)  and  invoking  Equation  (6)  yields 


N 


S(^ni  +  rlnj)  =  1  = 
n=l 


(12) 


thus, 


Substituting  Equation  (13)  into  Equations  (10)  and  (11)  gives 


(13) 


(14) 


(15) 


In  Equations  (14)  and  (15)  the  plus  sign  was  chosen  corresponding  to  the  direction  of  maximum 
function  increase.  By  changing  the  plus  to  a  minus  in  Equations  (14)  and  (15),  the  search  directions 
then  correspond  to  the  direction  of  maximum  function  decrease,  i.e.,  the  minus  sign  minimizes  the 
power  delivered  away  from  the  focus  or  tumor  site.  Partial  derivatives 


dFj  dFj 
dAnj'd*nj 


;n  =  1,2, •••,1V 


(16) 


represent  the  gradient  directions  for  maximum  function  increase. 

Because  F  is  measured  and  cannot  be  expressed  in  analytical  form,  the  partial  derivatives 
are  numerically  evaluated  using  finite  differences.  Thus 


dFj  _  &FAnj 
dA„j  2AAnj 


(17) 


dFj  _  AF*nj 
d*nj  2 


(18) 
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where,  as  shown  in  Figure  5, 


AF_ Anj  —  Fj(Anj  +  AAnj;$nj)  —  Fj(Anj  A  Anj ;  $nj )  ,  (19) 

AF$nj  =  Fj(Anj',<f?nj  +  A4*nj)  —  Fj(Anj;$nj  —  A$nj)  (20) 

and  AA„j  and  A 4>„j  are  the  maximum  step  sizes.  Assume  for  now  that  the  increments  AAnj  and 
A depend  on  the  iteration  number  j  and  transmit  element  index  n.  Substituting  Equations  (19) 
and  (20)  into  (17)  and  (18)  and  then  substituting  the  result  into  (14)  and  (15)  gives  the  desired 
result  for  the  search  directions 


Anj  — 


AFAn, 
A  An] 


rQnj  — 


AFytn, 
A  An] 


)2  +  ( 


AF*”> 

A$nj 


)21 


(21) 


(22) 


The  new  amplitude  and  phase  settings  of  the  (j  +  l)th  transmit- weight  configuration  are  computed 
according  to 


An,j+ 1  —  Anj  +  AAnjTyinj  (23) 

$n,j+ 1  =  $nj  4*  A$njTQnj  .  (24) 

For  the  current  software  implementation  of  the  gradient  search  in  these  experiments,  the  step  sizes 
are  assumed  (for  convenience)  to  be  independent  of  both  the  iteration  number  and  the  adaptive 
channel  number,  that  is 


AAnj  =  AA 

(25) 

A$nj  =  A$ 

(26) 

In  some  situations  it  may  be  desirable  to  change  the  step  size  at  each  iteration  as  considered 
in  Farina  and  Flam  [73]  and  Hasdorff  [72],  but  that  possibility  has  not  been  explored  in  current 
measurements. 
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Figure  5.  Figure  of  merit  with  transmit-weighi  dithering  for  optimum  search  directions. 
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2.3  System  Considerations 


Figure  6  is  a  block  diagram  of  an  adaptive  hyperthermia  system  controlled  by  the  gradient- 
search  algorithm.  7 .unsmit  weights  tojj,  •  ■  • ,  wn},-  •  ■ ,  w^j  at  the  jth  iteration  are  shown  at  the  top 
of  the  figure.  The  transmit  phased-array  antenna  induces  a  voltage  across  the  terminals  of  the  ith 
receive  field  probe  antenna.  For  any  given  configuration  of  the  transmit  weights,  each  weight  is 
dithered  by  a  small  amount  in  amplitude  and  phase,  and  the  received  powers  at  the  E-field  probes 
are  stored  in  the  computer  for  calculation  of  /,  search  directions,  and  updated  (j  +  l)th  transmit- 
weight  configuration.  One  transmit  weight  is  dithered  with  the  remaining  transmit  weights  in  their 
jth  state.  The  figure  of  merit  Fj  in  the  adaptive  hyperthermia  system  is  the  power  received  by  the 
auxiliary  probe  array,  as  indicated  in  the  block  diagram.  Search  directions  for  the  adaptive  transmit 
weights  are  based  on  maximizing  the  auxiliary  probe  array  received  power  and  are  computed  based 
on  Equations  (21)  and  (22).  Transmit  weights  for  the  next  configuration  (j  + 1)  are  computed  from 
Equations  (23)  and  (24).  The  adaptive  weight  vector  wa  is  achieved  when  the  ( j  +  l)th  weight 
configuration  has  converged.  Based  on  other  measurements  [10,12]  not  shown  here,  convergence 
occurs  typically  in  less  than  20  iterations,  depending  on  the  step  sizes  (A A  and  A4>). 
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ADAPTIVE  TRANSMIT  WEIGHT  ( J  +1)th  STATE 


Figure  6.  Gradient-search  algorithm  block  diagram  for  adaptive  focusing  hyperthermia 
system. 
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3.  MATERIALS  AND  METHODS 


The  hyperthermia  phased-array  system  used  in  these  measurements  is  the  Labthermics  Tech¬ 
nologies  Microtherm- 1000  planar  phased-array  applicator  (Figure  1),  using  a  15-  x  15-cm  planar 
array  with  a  4  x  4  grid  of  uniformly  spaced  square  waveguide  elements  operating  at  915  MHz 
(Figure  2).  Figure  7  is  a  photograph  of  the  planar  phased-array  applicator.  Due  to  the  water 
bolus  covering  the  applicator,  it  is  not  possible  to  see  the  actual  waveguide  radiators  of  the  ar¬ 
ray.  The  waveguide  aperture  dimensions  are  3.8  x  3.8  cm,  and  the  waveguides  are  loaded  with  a 
high-dielectric  constant  material  to  impedance-match  the  elements  at  915  MHz.  The  waveguides 
are  linearly  polarized  with  the  dominant  E-field  component  aligned  with  the  y-direction  (refer  to 
Figure  2).  The  variable  transmit  gain  and  phase  modules  are  controlled  by  digital- to-analog  (D/A) 
converters.  There  are  16  high-power  amplifiers  that  drive  the  waveguides  with  up  to  35-W  average 
power  per  channel.  The  power  applied  to  each  element  is  controllable  in  10%  increments  of  the  ad¬ 
justable  master  power  level.  Each  of  the  16  active  channels  has  an  electronically  controlled  variable 
phase  shifter  (8  bits,  0°  to  360°)  for  focusing  the  array.  The  phased-array  applicator  can  be  driven 
in  incoherent  mode  (with  16  independently  operating  sources  with  frequencies  near  915  MHz)  or 
coherent  mode  (with  a  single  frequency  generator  and  independent  control  of  the  relative  phase  of 
each  aperture  to  within  1.5°).  Invasive  Luxtron  [78]  fiber-optic  temperature  sensors  are  included 
with  the  Microtherm- 1000  system  to  monitor  the  temperatures  achieved  during  treatment. 

The  coupling  bolus  consists  of  deionized  water  contained  within  a  flexible  membrane,  which 
can  be  expanded  outward  to  6  to  8  cm  in  front  of  the  applicator  face  to  possibly  improve  the 
microwave  coupling  to  sharply  contoured  surfaces.  Bolus  temperature  is  controlled  by  circulating 
cooling  water  through  a  thermal  conduction  heat  exchanger  contained  in  the  outer  edges  of  the  bolus 
housing.  Cooling  capacity  and  transient  control  of  the  coupling  bolus  in  the  UCSF  Microtherm 
unit  were  improved  by  removing  the  plastic  tubing  heat  exchanger  inside  the  bolus  and  adding  a 
closed-circuit  pumping  system  to  circulate  temperature-controlled  deionized  water  directly  through 
the  bolus  compartment  in  series  with  a  heat  exchanger  mounted  in  a  temperature-regulated  water 
bath.  A  cool  water  bolus  between  the  patient  and  the  phased  array  prevents  excess  heating  of  the 
skin  surface.  The  bolus  is  filled  with  circulating  deionized  water,  which  has  a  very  low  microwave 
propagation  loss. 

The  UCSF  hyperthermia  array  (Figure  8)  is  controlled  and  monitored  by  an  MS-DOS-based 
personal  computer  system,  which  implements  the  gradient-search  adaptive  focusing  algorithm. 

The  E-field  probe  [59,60,62]  used  in  these  measurements  was  fabricated  at  UCSF  and  inte¬ 
grated  with  the  Microtherm- 1000  system.  The  E-field  sensor  is  a  1.9-mm  OD  Schottky  detector 
diode  (HP-3486),  operating  in  the  square  law  region  with  the  diode  leads  arranged  to  form  a  dipole 
antenna  of  1-cm  length.  The  sensor  is  coupled  to  an  amplifier  circuit  (Figure  9)  via  high-resistance 
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Figure  7.  The  915-MHz  microwave  planar  phased-array  applicator  used  in  the  Labther- 
mics  Microtherm- 1000  system. 


lead  material  (60  kfl/m)2  to  minimize  perturbation  of  the  fields.  The  rectified  DC  output  voltage  is 
a  function  of  the  E-field  squared  in  the  orientation  of  the  dipole  leads.  The  dipole  piobe  hardware 
and  probe  (x ,y,z)  positioner  are  shown  in  Figure  10.  A  12-bit  analog-to-digital  (A/D)  converter 
samples  the  E-field  probe  signal. 

A  liquid-type  (salt  and  sugar  dissolved  in  deionized  water)  muscle-equivalent  phantom  was 
used  ( a  —  1.38  S/m,  relative  permittivity  er  =  54.7  [79]),  having  electrical  properties  similar  to 
high  water  content  tissues  at  915  MHz  (cr  =  51.0  and  a  =  1.28  S/m  [80]).  The  liquid  phantom  was 
placed  within  a  50-  x  50-  x  22-cm-deep  Plexiglas  tank  (Figure  10). 


2From  Holaday  Industries;  1.2-mm  OD,  including  outer  insulation. 
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Figure  8.  UCSF/MIT  adaptive  hyperthermia  experiments  equipment. 


JIT51M 

HIGH-RESISTANCE  LEADS 


Figure  9.  Amplifier  circuit  for  E-field  probe. 
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A  computer  program  that  implements  the  gradient-search  algorithm  was  developed  by  MIT 
Lincoln  Laboratory  and  integrated  by  UCSF  with  the  Labthermics  system.  The  MIT/UCSF  in¬ 
tegrated  software  allows  the  operator  to  run  an  adaptive  focusing  algorithm  that  uses  the  output 
power  of  the  E-field  probe  as  a  feedback  signal  during  the  gradient  search.  The  adaptive  focusing 
algorithm  uses  only  phase  control  to  maximize  the  focal  point  E-field. 


2079M-1 


Figure  10.  Dipole  probe,  probe  (x,y,z)  positioner  system,  and  the  muscle-equivalent 
phantom  tank. 
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4.  MEASURED  RESULTS 


4.1  Review  of  Earlier  Investigations  at  UCSF 

This  section  reviews  data  previously  measured  by  UCSF  that  demonstrate  several  capabilities 
of  the  Labthermics  system.  At  UCSF,  initial  efforts  to  study  the  electrical  phasing  capabilities  of  the 
Microtherm- 1000  unit  proved  difficult  due  to  irregular  radiation  beam  profiles  and  magnitudes  of 
individual  waveguide  apertures  [81]  and  the  inability  to  accurately  calibrate  or  preset  the  individual 
phase  shifters  to  absolute  values.  Attempts  at  using  theoretically  selected  phases  that  were  applied 
to  each  element  were  unsuccessful  or  produced  marginal  results  in  a  water  phantom.  To  compensate 
for  these  problems  and  tc  accommodate  expected  tissue  heterogeneities  in  vivo,  a  phase  selection 
or  optimization  scheme  was  developed  to  iteratively  select  the  phasing  of  each  transmit  channel  to 
maximize  the  measured  SAR  at  the  desired  location.  A  computer/software  system  linked  to  the 
Microtherm  unit  turns  on  individual  elements  of  the  phased  array  in  a  predetermined  sequence  and 
iteratively  adjusts  the  phasing  until  a,  maximum  SAR  is  recorded.  The  next  element  in  the  sequence 
is  turned  on  and  phase-adjusted  similarly  to  maximize  the  SAR.  This  phasing  process  continues 
until  all  elements  are  phase-focused.  This  phasing  approach  attempted  to  focus  in  deionized  water, 
tissue-equivalent  liquid  phantoms,  and  in  vivo  pig  thigh  [58].  With  the  UCSF  phasing  scheme,  it 
was  possible  to  focus  the  Microtherm  at  depths  of  7  cm  in  the  deionized  water  phantom  as  shown 
in  Figure  11(a).  (The  measured  E-field  radiation  pattern  contours  in  Figure  11  are  given  in  10% 
intervals.)  Clearly,  in  the  deionized  water  maximum  radiation  occurs  in  the  vicinity  of  the  desired 
focus.  Based  on  the  50%  SAR  value,  the  maximum  penetration  depth  is  15  cm  in  the  deionized 
water  medium.  Using  the  same  algorithm,  similar  focusing  attempts  in  muscle-equivalent  phantom 
[Figure  11(b)]  did  not  generate  a  useful  focused  beam  but  produced,  in  essence,  a  collimated  beam 
with  a  beam  peak  effectively  at  the  surface  of  the  phantom.  The  data  in  Figure  11(b)  indicate  a 
maximum  penetration  depth  of  3  cm,  based  on  the  50%  contour  that  reaches  z  =  3  cm.  (Microwave 
propagation  loss  in  deionized  water  is  0.3  dB/cm,  whereas  in  muscle-equivalent  phantom,  loss  is 
close  to  3  dB/cm.)  Propagation  loss  in  muscle  tissue  clearly  contributes  greatly  to  the  inability  to 
focus  the  915-MHz  planar  phased  array  at  the  desired  depth.  Even  with  the  degradation  in  focus, 
a  3-cm  penetration  depth  is  potentially  useful  for  certain  tumors. 

The  efficacy  of  the  UCSF  focusing  algorithm  was  evaluated  in  a  pig  thigh  model,  where  the 
phase  was  optimized  using  an  implanted  SAR  sensor.  Figure  12  shows  the  in  vivo  experimental 
set-up  with  trial  focus  positions  (z  =  2.5  and  4.5  cm)  indicated  by  the  shaded  circles  located  on  the 
central  axis.  In  the  incoherent  mode  with  the  16  waveguides  operating  independently,  the  measured 
radiation  pattern  in  Figure  13(a)  indicates  a  maximum  penetration  depth  of  about  1  cm  on  the 
central  axis,  consistent  with  a  microwave  propagation  loss  of  3  dB/cm.  Figures  13(b)  and  (c)  show 
the  measured  radiation  pattern  for  a  focus  at  z  —  2.5  and  4.5  cm,  respectively.  Focused  at  z  —  2.5 
cm  [Figure  13(b)],  the  measured  maximum  penetration  depth  (on  the  principal  axis)  is  2.3  cm,  and 
the  beam  peak  is  at  z  =  1  cm.  Focused  at  z  =  4.5  cm  [Figure  13(c)],  however,  penetration  depth 
is  only  1.9  cm,  and  the  beam  peak  is  at  the  surface  of  the  target.  Intuitively,  one  would  expect  the 
z  =  4.5-cm  focus  to  penetrate  deeper  than  the  z  =  2.5-cm  focus.  Instead,  measurements  indicate 
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UCSF  MEASUREMENTS  im<M 
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Figure  11.  Measured  E-field  radiation  pattern  for  Microtherm- 1000  planar  phased  array 
focused  at  (a)  z  —  7  cm  m  deionized  water  phantom  (propagation  loss  =  O.S  dB/cm), 
resulting  in  a  well-focused  beam;  (b)  z  =  4  cm  in  muscle-equivalent  liquid  phantom  (prop¬ 
agation  loss  —  S  dB/cm),  resulting  in  a  beam  that  is  not  centered  at  the  desired  focusing 
depth. 


22 


a  reduction  in  penetration  depth  when  the  focus  is  increased  beyond  z  =  2.5  cm.  This  discrepancy 
was  the  principal  reason  that  UCSF  was  interested  in  comparing  the  MIT  Lincoln  Laboratory 
gradient-search  algorithm  with  the  systematic  iterative- search  algorithm  currently  used  with  the 
Microtherm  system.  Figure  14  compares  the  measured  SAR3  and  temperature  distributions  for  the 
z  —  2.5-cm  focus  case,  and  the  two  sets  of  data  are  consistent;  that  is,  the  peak  temperature  value 
occurs  near  the  peak  of  the  SAR  distribution. 
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Figure  12.  UCSF  experimental  test  set-up  for  «n  vivo  measurements  of  pig  thigh. 


From  these  measurements  it  is  concluded  that  by  optimizing  the  phasing  of  the  Microtherm  it 
is  possible  to  focus  in  a  low-loss  deionized  water  phantom.  The  focused  power  depositions  in  muscle 
phantom  and  tissue  do  not  produce  a  useful  focus  but  illustrate  a  possible  improvement  in  pene¬ 
tration  depth.  One  question  that  remains  concerns  the  validity  of  the  UCSF  optimization  routine 


3Defined  here  as  SAR  =  cAT/At  where  AT  is  the  change  in  temperature  over  time  interval  At  and 
c  is  the  specific  heat. 
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UCSF  MEASUREMENTS 


•mi* 4 


Figure  13.  Measured  E-field  radiation  pattern  for  Microtherm- 1000  planar  phased  array 
illuminating  muscle-equivalent  liquid  phantom  (propagation  loss  =  3  dB/cm):  (a)  inco¬ 
herent  mode  (all  elements  incoherently  driven);  focused  at  z  =  (b)  2.5  and  (c)  f.5  cm. 
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and  whether  better  focusing  in  a  lossy  medium  could  be  obtained  with  an  adaptive  optimization 
algorithm  as  considered  previously  [7,10,12].  New  measurements  follow. 

4.2  New  Measurements  at  UCSF 

UCSF  and  MIT  Lincoln  Laboratory  began  collaborating  in  June  1992  to  implement  an 
adaptive  focusing,  gradient-search  algorithm  with  the  Microtherm- 1000.  Experiments  using  the 
UCSF  and  MIT  iterative-search  algorithms  were  performed  in  the  deionized  water  and  the  muscle- 
equivalent  liquid  phantoms.  Phase  focusing  was  attempted  in  the  central  axis  of  the  Microtherm 
applicator  at  depths  of  6  and  8  cm  from  the  lowest  ridge  on  the  applicator  housing  with  an  ad¬ 
ditional  3-cm  path  length  in  deionized  water  from  the  apertures  to  this  point.  Measured  E-field 
patterns  (Figure  15)  with  the  focus  set  at  8  cm  indicate  that  both  optimization  schemes  work  simi¬ 
larly,  producing  well-defined  focused  beams  at  the  desired  depths.  Current  measurements  have  not 
considered  field  mapping  in  other  vertical  and  horizontal  planes  that  are  necessary  to  characterize 
the  locations  and  magnitudes  of  side-  and  possible  grating  lobes.  It  is  interesting  to  note  that 
the  radiation  patterns  in  Figure  15  exhibit  a  z-dependent  ripple  with  spacing  equal  to  «  1.8  cm. 
This  ripple  is  likely  due  to  a  standing-wave  pattern  created  by  incident  and  reflected  waves  in  the 
phantom  water  tank.  Assuming  a  dielectric  constant  of  80  and  an  electrical  conductivity  of  0.19 
S/m,  the  wavelength  of  915-MHz  microwave  radiation  in  deionized  water  is  3.6  cm.  Thus  the  ripple 
spacing  is  equal  to  one- half  wavelength  as  expected  with  a  standing- wave  pattern. 

Measurements  in  the  muscle-equivalent  phantom  were  performed  at  a  depth  of  4  cm  for 
three  cases  of  amplitude  illumination:  equal  applied  at  100%  of  selected  power  level  (Figure  16), 
adjusted  uniform  using  premeasured  SAR  amplitudes  for  each  element  (Figure  17)  [81],  and  inverse- 
tapered  by  applying  double  the  power  to  the  outer  elements  (Figure  18).  These  radiation  pattern 
measurements  were  similar  to  those  obtained  in  earlier  published  studies  [58].  Notice  that  there 
are  only  minor  differences  between  the  three  radiation  patterns  of  Figures  16  to  18.  It  is  likely 
that  the  four-element  group  at  the  center  of  the  array  is  the  dominant  contributor  to  the  radiation 
pattern  shape.  Shown  last  is  the  convergence  of  the  method  of  steepest- ascent  algorithm  for  the 
case  of  the  4-cm  focal  depth.  Appendix  A  lists  the  raw  A/D  data  for  the  measured  E-field  probe 
power  in  an  experiment  consisting  of  15  iterations  of  the  phase  focusing  steepest- ascent  algorithm. 
The  A/D  converted  data  are  changed  to  power  in  decibels  by  computing  101og(p),  where  p  is  the 
value  (measured  power)  of  the  A/D  converter.  The  measured  E-field  at  the  focus  as  a  function 
of  iteration  number  is  displayed  in  Figure  19.  Focal-probe  power  increases  by  17  dB  over  the  15 
iterations  with  the  data  indicating  convergence  in  about  7.  [Note:  The  maximum  step  size  given 
by  Equation  (26)  is  five  D/A  states  at  each  iteration.  Also,  each  iteration  takes  less  than  1  min  to 
execute,  thus  the  array  is  automatically  focused  in  less  than  10  min.] 

The  results  from  these  investigations  indicate  that  the  low  number  (16)  of  array  antenna 
elements,  high  attenuation  of  the  signal  from  outer  elements,  and  irregular  beam  patterns  of  the 
individual  elements  dictate  that  the  Microtherm  unit  cannot  produce  a  well-defined  focus  at  any 
appreciable  depth  beyond  3  cm  in  lossy  muscle  tissue.  Beam  shaping  and  possible  improvement  in 


25 


Figure  H .  Comparison  of  (a)  measured  SAR  and  (b)  temperature  distributions  (steady 
state,  10  min)  for  Microtherm- 1000  system  focused  at  2.5  cm  in  pig  thigh. 
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Figure  15.  Comparison  of  measured  two-dimensional  radiation  patterns  of  the 
Microtherm- 1000  with  focal  depth  of  8  cm  in  deionized  water  phantom:  (a)  UCSF  system¬ 
atic  iterative-search  focusing  algorithm  and  ( b )  MIT  Lincoln  Laboratory  adaptive  gradient- 
search  focusing  algorithm. 
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MIT  GRADIENT  SEARCH.  MUSCLE  PHANTOM  (29  JUNE  1992) 
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Figure  16.  Measured  two-dimensional  radiation  pattern  of  the  Microtherm-1000  using 
MIT  Lincoln  Laboratory  adaptive  focusing  algorithm  with  focal  depth  of  4  cm  in  muscle- 
equtvaleni  liquid  phantom.  Power  level  of  each  waveguide  is  at  100%  relative  to  each  other 
(uniform  illumination). 


penetration  depth  with  different  illumination  strategies  is  one  aspect  that  can  be  studied,  but  sig¬ 
nificant  improvement  is  not  expected.  Curved  phased-array  applicators  with  the  radiating  antenna 
elements  pointed  toward  the  focus  may  improve  effective  penetration  depth  and  will  be  investigated 
in  future  measurements.  Because  two  proven  types  of  phase  optimization  routines  yield  similar  re¬ 
sults,  confidence  exists  that  the  measured  data  truly  depict  the  best  possible  performance  of  the 
Microtherm  system. 
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MIT  GRADIENT  SEARCH,  MUSCLE  PHANTOM  (29  JUNE  1992) 
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Figure  17.  Measured  iwo-dimensional  radiation  pattern  of  the  M icro therm- 1 000  using 
MIT  Lincoln  Laboratory  adaptive  focusing  algorithm  with  focal  depth  of  4  cm  in  muscle- 
equivalent  liquid  phantom.  Power  level  of  each  waveguide  is  adjusted  to  provide  equal 
amplitude  as  measured  at  the  focusing  probe  (adjusted  uniform  illumination). 


29 


zimt-io 


50 

60 

20 

60 

100 

40 

y 

40 

100 

100 

50 

30 

80 

60 

60 

30 

80 

ARRAY 

ILLUMINATION  (%) 
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Figure  18.  Measured  two-dimensional  radiation  pattern  of  the  Microtherm- 1000  using 
MIT  Lincoln  Laboratory  adaptive  focusing  algorithm  with  focal  depth  of  I  cm  in  muscle- 
equivalent  liquid  phantom.  Power  level  of  the  outer  ring  of  waveguides  is  adjusted  to 
provide  twice  the  power  of  the  inner  ring  of  elements  as  measured  at  the  focusing  probe 
(inverse  taper  illumination ). 
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UCSF/MfT  MEASUREMENTS  (29  JUNE  1992) 


ITERATION  NUMBER 


Figure  19.  Measured  E-field  focused  at  z  =  4  cm  as  a  function  of  adaptive  phase  focusing 
gradient-search  iteration  number  for  a  muscle- equivalent  liquid  phantom.  Power  level  at 
the  focus  increases  by  17  dB  due  to  the  adaptive  gradient  search. 
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5.  CONCLUSION 


An  adaptive  focusing,  915-MHz,  microwave  planar  phased-array  hyperthermia  system  has 
been  experimentally  demonstrated  by  modifying  a  commercial  version  at  the  University  of  California 
at  San  Francisco.  The  measurements  shown  in  this  report  demonstrate  that  the  E-field  radiation 
pattern  of  a  hyperthermia  planar  phased  array  can  be  controlled  by  an  adaptive  feedback  gradient- 
search  focusing  algorithm.  E-field  focusing  (power  maximization)  at  a  desired  phantom  target 
position  has  been  experimentally  demonstrated. 

In  this  report,  the  adaptive  algorithm  is  a  gradient-search  feedback  technique  (method  of 
steepest  ascent)  derived  in  Section  2.  The  gradient  search  employs  transmit- weight  phase  dithering 
with  E-field  probe  power  measurements  at  the  desired  target  position  to  calculate  the  required 
gradient-search  directions  for  microwave  focusing.  The  steepest  ascent  algorithm  produces  two- 
dimensional  radiation  patterns  that  are  similar  to  those  produced  by  a  systematic  iterative  focusing 
algorithm  previously  demonstrated  with  measurements  at  UCSF. 

E-field  focusing  is  intended  to  maximize  the  microwave  power  delivered  to  a  tumor  site.  The 
data  presented  in  this  report  suggest  that  this  goal  may  be  achieved  for  shallow  tumors  (less  than 
3  cm  beneath  the  skin  surface)  with  an  adaptive  focusing,  hyperthermia  microwave  planar  phased- 
array  system.  Measurements  have  proven  that  a  planar  phased  array  at  915  MHz  can  only  provide 
a  useful  focus  in  lossy  (3  dB/cm)  muscle  tissue  at  a  maximum  depth  of  3  cm. 

Because  of  the  hard-  and  software  modifications  made  at  UCSF,  the  16- channel  915-MHz 
Microtherm- 1000  system  can  serve  as  a  test  bed  for  new  types  of  adaptive  phased -a, -ray  applicators. 
For  example,  it  may  be  possible  to  generate  a  focus  deeper  than  3  cm  with  a  curved  phased-array 
geometry  because  radiating  antennas  on  a  curved  contour  tend  to  have  peak  radiation  toward 
the  focus  and  propagation  path  length  between  the  antennas  and  the  focus  will  vary  slightly  from 
element  to  element.  Thus  the  curved  (compared  with  the  planar)  array  may  provide  a  better  focus  at 
depth.  Future  measurements  will  investigate  a  new  noninvasive,  monopole,  adaptive  phased-array, 
915-MHz  microwave  applicator  (recently  designed  and  fabricated  at  MIT  Lincoln  Laboratory)  with 
other  types  as  they  become  available. 
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APPENDIX  A 

Measured  E-F:eld  Probe  Raw  Data 


Measured  Raw  A/D  Converter  Data  for  E-Field  Probe  at  Focus 
Adaptive  Focus,  Muscle-Equivalent  Liquid  Phantom 

Focus  at  z  *  4  cm  of  Phantom  (run  wfa8z).  Method  of  Steepest-Ascent  Algorithm 

j  *  iteration  number 
PI  *  probe  1  (tumor  site) 

j  PI 


0  1.2 

1  4.6 

2  10.1 

3  18.5 

4  27.0 

5  35.4 

6  42.6 

7  48.5 

8  52.8 

9  55.8 

10  57.9 

11  58.8 

12  61.1 

13  62.2 

14  62.8 

15  64.2 

Note:  The  above  A/D  measured  data  are  converted  to  decibels  by  taking 
10  times  the  log  of  the  tabulated  numbers . 
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